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1. Introduction 



Over the last two decades theory and experiment on quantum disordered and 
chaotic systems have been an extremely successful field of research. The mile- 
stones on this rout were: 

• Weak Anderson localization in disordered metals 

• Universal conductance fluctuations 

• Application of random matrix theory to quantum disordered and chaotic sys- 
tems 

However the mainstream of research was so far limited to a linear response of 
quantum systems where the kinetic and thermodynamic properties are calculated 
essentially at equilibrium. The main goal of this course of lectures is to develop a 
theory of nonlinear response, to a time-dependent perturbation in the same way as 
the theory of Anderson localization and mesoscopic phenomena. It will require 
an extension of an existing methods to non-equilibrium phenomena. 
The four lectures will include the following topics: 

• Perturbative theory of weak Anderson localization 

• Keldysh formulation of nonlinear response theory 

• Mesoscopic rings under AC pumping and quantum rectification 

• Theory of weak dynamic localization in quantum dots 



2. Weak Anderson localization in disordered systems 

In this section we consider the main steps in the calculus of weak localization 
theory [1,2]. The main object to study is the frequency-dependent conductivity 



where Gj T are retarded (advanced) electron Green's functions, /(e) is the 
Fermi energy distribution function and v a is the velocity operator . One can 




(2.1) 
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convince oneself using the representation in terms of exact eigenfunctions ^> n (r) 
and exact eigenvalues E n 

r R/A _^ W„(r)W;(rQ 

G <'^-^e-E n ± l (22) 



that Ea. J2.1> reduces to a familiar Fermi Golden rule expression: 

o a pW) = 2tt ^ (n|J a |m)(m| Jp\n) 5(E n - E m - u) ff" ( 2 -3) 

B„>£ m ™ m 



Ea. (l2.2> is convenient to prove exact identities but is useless for practical calcu- 
lus. In disordered systems with the momentum relaxation time r the following 
expression for the disorder average Green's function is very useful: 

(G R/ \ = _ f L , J _ - (2-4) 

where £ p = e(p) — Ep is related with the electron dispersion law e(p) relative 
the Fremi energy Ef. In what follows we will assume £ = £ p depending only 
on |p|. Then £ and the unit vector of the direction of electron momentum n 
constitute a convenient variables of integration over p: 

.dp -> J J da... (2.5) 

For a typical metal with Ef far from the band edges and the energy scale of 
interest e << Ep the density of states v(£) can be approximated by a constant 
v = i/(0) at the Fermi level and the integration over £ can be extended over the 
entire real axis (— oo, +oo). 



2.1. Drude approximation 

The simplest diagrams for the frequency-dependent conductivity are shown in 
Fig.l. The diagrams Fig. la,b are given by the integrals 

f f + °° d£ 
dnv a v/3 — 1 — — . - = (2.6) 

where v a ~ VFn a . 

It is important that the integral over £ is the pole integral with all poles in 
the same complex half-plane. This is why such integrals are equal to zero. In 
contrast to that the diagrams of Fig.lc,d that contain both retarded (G R ) and ad- 
vanced (G A ) Green's functions correspond to a similar integral with poles lying 
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(a) (b) (c) (d) 



Fig. 1 . Diagrams for the Drude conductivity. 



in different half-planes of £. Such an integral over £ can be done immediately 
using the residue theorem: 



« ^ (2.7) 



The angular integral is trivial: 



/ 



dnv a vp = v% (2.8) 



where d is the dimensionality of space. 

Now all what we need to compute a a p using the simplest diagrams of Fig.l 
is the identity: 



/ 



+ 00 

[f(e)-f(e-u)}de = -Lu. (2.9) 



This identity holds for all functions /(e) obeying the Fermi boundary conditions 
f(e) —* at e — > +oo and /(e) — > 1 at e — > — oo. 

The result of the calculations is the Drude conductivity: 

(D) e 2 vD 

where Z? = vj^r/d is the diffusion coefficient. 
2.2. Beyond Drude approximation. 

The diagrams that were not included in the Drude approximation of Fig.l con- 
tain cross-correlations between the exact electron Green's functions G R and G A 
which arise because both of them see the same random impurity potential U (r). 



8 



We make the simplest approximation of this potential as being the random Gaus- 
sian field with zero mean value and zero-range correlation function: 

6(r - r'l 

(U(r)U(r')) = -L 1, (2.11) 

LHVT 

The simplest diagrams beyond the Drude approximation are the ladder series 
shown in Fig. 2. The dotted line in these diagrams represent the momentum 
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Fig. 2. Ladder diagrams. 



Fourier transform of correlation function (U(r)U(r')) which is aconstant 1/2-kvt 
independent of the momentum transfer p — p'. This makes the quantities v a and 
v'p completely independent of each other. Then the angular integration 

J dnv a =Q (2.12) 

results in vanishing of all the ladder diagrams of Fig. 2. 

The next diagram with cross correlations contain an intersection of dotted 
lines (see Fig. 3). 

In contrast to a diagram with two parallel dotted lines where all integrations 
over momenta are independent, here the momentum conservation imposes a con- 
straint: 

P + p' = Pi+P2- (2.13) 

At the same time the form Ea. (l2.4> of the average Green's function suggests that 
the main contribution to the momentum integration is given by momenta confined 
inside a ring of the radius p p and the width 1 /£ <C pf where I is the elastic mean 
free path £ — vpr [see Fig.4a]. The constraint Ea. J2.13t implies that not only p 
but also p + Pi + P2 = p' should be inside a narrow ring [see Fig. 4b]. Thus 
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Fig. 3. The first of the "fan" diagrams. 




Fig. 4. Regions that make the main contribution to momentum integrals. 



the effective region of p integration is an intersection of two narrow rings which 
volume is reduced by a large factor of {pf() compared to an unconstrained case. 

2.3. Weak localization correction. 

The lesson one learns from the above example is that any intersection of dotted 
lines in systems with dimensionality greater than one brings about a small pa- 
rameter l/(pp£). However this does not mean that all such diagrams should be 
neglected. As a matter of fact they contain an important new physics of quantum 
coherence that changes completely the behavior of low dimensional systems with 
d = 1, 2 at small enough frequencies ui leading to the phenomenon of Anderson 
localization. 

For this to happen there should be something that compensates for the small 



10 



factor l/(pjr£). We will see that this is a large return probability of a random 
walker or a particle randomly scattered off impurities in low-dimensional sys- 
tems. On the formal level the corresponding diagrams are just the multiple scat- 
tering "fan" diagrams shown in Fig. 5. One can see using the momentum con- 
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Fig. 5. Infinite series of "fan" diagrams. 



servation that such diagrams have only one constraint of the type Eq.( I2.13> in- 
dependently of the number of dotted lines. This statement becomes especially 
clear if one rewrites the fan series in a form Fig. 5b which contains a ladder series 
with all dotted lines parallel to each other and only one intersection of solid lines 
representing electron Green's functions. It is this ladder series called Cooperon 
which describes quantum interference that leads to Anderson localization. 

To make connection between random walks (or diffusion) in space and a quan- 
tum correction to conductivity and to see how a compensation of the small pa- 
rameter l/(pp£) occurs we compute the Cooperon C(k, u>) for a small value of 
Pi + P2 = k [see Fig. 6.]. This series is nothing but a geometric progression with 
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Fig. 6. Summation of the ladder series for a Cooperon. 



the first term 

Ho - j d V G*(j>)Gf_ u {- V + K) (2.14) 
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and the denominator q = 2irvTHo. Thus we have 

C(k,w) = — — . (2.15) 

i - q 

To compute Ilo we write 

Ho = , n - v> / ^ / + — — — . (2.16) 



(2tti/t) 2 J ( e - -u-t + vpnk- 

In Ea.( l2.16> we used an approximation £- p +k ~ £p — w^nk which is valid as 
long as |k| <C pf- The pole integral in Ea. (l2.16> can be done immediately with 
the result: 

q= [ da- . * -. , . (2.17) 

J 1 — IUJT + IIXIK 

One can see a remarkable property of Ea. (l2.17> : the denominator of the geomet- 
ric series Ea. J2.15> is equal to 1 in the limit ujt, \k\£ — ► 0. At finite but small 
cut, \k\£ <C 1 one can expand the denominator of Eq. j2.^ to obtain: 

q ^l + iuJ T-£ 2 k 2 /d; C(k,u) « — ^ 1 . . (2.18) 

We immediately recognize an inverse diffusion operator in C (k, cj) which is di- 
vergent at small u> and k. This divergence is the cause of all the peculiar quantum- 
coherence phenomena in systems of dimensionality d < 2. 

In particular the quantum correction to conductivity given by the diagram of 
Fig. 5 can be written as: 

^ = ^x^£^C(k, W ). (2.19) 

k 

where T a @ is the "Hikami box" shown in Fig. 7: Its analytic expression is given 
by: 

J dpv Q (p)v (-p+k) Gf(p)G^_ w (p)Gf(-p+k)G^(-p+k).(2.20) 

In the limit ujt <C 1 and \k\£ <C 1 one can set u> = k = in Eq.( I2.20> . Then 
v(— p + k) = — v(p) = -vpii and after the pole integration over £ and the 
angular integration over n we obtain for r a ^: 



/ +00 ( £ - C+ J^-g-^) 2 = " 4m/T2j °- (2 ' 21) 
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V A 




(a) 



(b) 



Fig. 7. The diagram representation for the Hikami box (a) and for the quantum correction to conduc- 
tivity (b). 

Finally using Eq.( l2.19> we get: 



This is the celebrated formula for weak Anderson localization. One can see that 
the quantum correction is given by the sum over momenta of the diffusion prop- 
agator, that is, it is proportional to the return probability at a time ~ 1/u) for a 
random walker in the d-dimensional space. A remarkable property of random 
walks in low-dimensional space is that the return probability increases with time. 
That is why the quantum correction to conductivity increases with decreasing 
the frequency ui as in a quasi-one dimensional wire and as log(l/o>) in a 

two-dimensional disordered metal. 

The structure of Ea. d2.22> suggests a qualitative picture of weak Anderson 
localization. It is an interference of two trajectories with a loop that differ only in 
the direction of traversing the loop [Fig.8]. Although the phase that corresponds 
to each trajectory is large and random, the phase difference between them is zero 
because of the time-reversal invariance. As a result they can interfere and amount 
to a creation of a random standing wave which is the paradigm of Anderson 
localization. 



3. Non-linear response to a time-dependent perturbation 

In this section we review the main steps of the Keldysh formalism [3,4] which 
are necessary to describe a quantum system of non-interacting electrons subject 
to external time-dependent perturbation V(t). 



\ TXV Vol ^— ' 

\ k 




(2.22) 
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3.1. General structure of nonlinear response function. 
The matrix Keldysh Green's function 



G=( n ^1 (3-D 



G R G K 
G 

contains - besides familiar retarded and advanced Green's functions G 11 = G R 
and G 22 = G A - also the third, Keldysh function G 12 = G K . The latter is 
the only one needed to compute an expectation value of any operator O both in 
equilibrium and beyond: 

O = -iTr(OG K ). (3.2) 

The retarded and advanced Green's functions appear only at the intermediate 
stage to make it possible to write down Dyson equations in the matrix form. It is 
of principal importance that the component G 21 is zero. One can show that this 
is a consequence of causality [5]. In equilibrium there is a relationship between 
G K and G R I A . It reads: 



G 



K 



(G R -G A ) E t & nh(^L). (3.3) 
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As without time-dependent perturbation the system is in equilibrium, the corre- 
sponding components Gq^ A and Gq e of the matrix Keldysh Green's function 
obey the relationship Ea. (l3.3> . 

In the presence of the perturbation one can write the Dyson equation: 

G = G + G VG, (3.4) 

where the operator of external time-dependent perturbation V(£) is proportional 
to a unit matrix in the Keldysh space. The structure of perturbation series for G K 



V V V V 




V V V V 



R-A junction 



V [tanh( e ,/2T) -tanhfr 10 )\ 



Fig. 9. The diagram representation of the "anomalous" term in the Keldysh function. 



that corresponds to Ea. (l3.4> is as follows 



GK y^ll 
' ' i 



-G^VGtf 



-Gl 2 VG 2 2 



-G^V ...G^VGl'VG 22 



■G n 



.(3.5) 



As Gf = and G\ 2 = (G§ - G5 4 ) tanh(e/2T) each term of the perturba- 
tion series Eq. (I3.5> is a string of successive Gq functions followed by a string 
of Gq functions with only one switching point between them where the factor 
tanh(e/2T) is attached to. For the same reason the components G R and G A are 
the series that contain only Gq or Gq, respectively: 



G R/A = Q R/A + G R/Ay G L 



R/A 



Gq /a VGq /a VG 

R/A 



R/A 



(3.6) 



Using Ea. d3.6t one can replace the string of G functions in Ea. P.5> by the 

R/A 

exact retarded or advanced Green's function G £ ' , that depends on two energy 
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variables because of the breaking of translational invariance in the time domain 
by a time dependence of perturbation. The result is: 

G£, = G R e , tanh (j^ - tanh (^) G^, + (3.7) 

The first two terms in Eq. j3^ just reproduce the structure of the equilibrium 
Keldysh function Eq.( I3.3> . The most important for us will be the last, so called 
"anomalous" term G anonl in Ea. J3.7> as it contains both G R and G A functions 
which makes nonzero the pole integrals over £ and allows to build a Cooperon. 
This term is graphically represented in Fig. 9. 

One can see that of all the vertices with a perturbation operator V one is spe- 
cial: it is a switching point between the strings of retarded and advanced Green's 
functions, so called R — ^junction [5]. It is convenient to switch from the energy 
to the time domain. In this representation the anomalous part of G K reads [5]: 

G%T = Gf.M Gtf (V(h) - V{ti)) &o(*i - *a). (3-8) 

In Eq. (I3.8> we denoted the Fourier transform of tanh(e/2T) as: 

f dp ?T 
ho® = / 7T e tEt tanh( £ /2T) = . (3.9) 
J ztt smn(7ri t) 

We also assume integration over repeated time variables. 

The graphic representation of Eq. d3~8l is given in Fig. 10. Eqs. j3.7H3.8> give 
a general structure of nonlinear response to a time-dependent perturbation. The 
main feature is the distinction between a regular V vertex that connects two 
Green's functions with the same analytical properties (both G R or both G A ) and 
the R — A junction. One can see that the time variables do not match at the R- 
A junction because of the attached energy distribution function, while in at any 
regular vertex V they do match. 

3.2. Approximation of single photon absorption/emission 

While an exact solution of the nonlinear response is a difficult task, there is an 
essentially nonlinear regime where a regular solution can be found. The corre- 
sponding approximation consists of expanding each disorder averaged retarded 
or advanced Green's function up to the second order in V: 

(G R / A ) « (G R/A ) + (G R/A )V(G R/A ) + (G R/A )V(G R/A )V(G R/A )(3A0) 



tanh 



f- 

V2T 



tanh 



El — UJ 
2T 
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G R (t,tD 



ti 

t 2 



G A (t 2 ,t') 



h(t 1 -t 2 )[V(t 2 )-V(t 1 )] 
non-locality in time 



Fig. 10. Anomalous contribution to G K in the time domain. 



Yet, since the ladder series (a Cooperon) contains an infinite number of retarded 
(advanced) Green's functions, this approximation leads to an essentially non- 
linear results that do not reduce to a finite order response function. We will 
see that this approximation implies a sequential absorption/emission of photons 
rather than the multiple -photon processes. 

Another approximation which is similar to lot <C 1 is based on the fact that 
the disorder average Green's function Eq.( I2.4> decays exponentially in the time 
domain: 

(Go /A ) f =TiO(±t)e~ u ^e- t/T . (3.11) 

Assuming that the momentum relaxation time is much shorter than the character- 
istic scale (e.g. the period) of a time-dependent perturbation one can approximate 
(G Q )t by a <5-function and its derivative: 

(< M ) t « S(t) (< M ) £=0 - id t 5(t)(G« /A ) 2 £=0 . (3.12) 



4. Quantum rectification by a mesoscopic ring. 

To proceed further on we need to make some assumptions about the time-dependent 
perturbation and specify an observable of interest. 

In this section we consider a problem of rectification of an ac signal by a 
disordered metal. It has first been considered in Ref. [9] in a single-connected 
geometry where the entire effect is due to mesoscopic fluctuations. Here we 
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focuse on the case of a quasi-one dimensional disordered metal ring pierced by 
a magnetic flux <fi(t) that contains both a constant part <fi and an oscillating part 
4>ac(t) [6-8]. In this case the topology of a ring and the presence of a constant 
magnetic flux makes it meaningful to study the disorder-average rectification 
effect. 

The time-dependent perturbation in this case is: 

V = -vMt), ¥»(t) = ^%^, (4.1) 

where L is the circumference of the ring, and (f>o — he/ e is the flux quantum. 
We assume that the ring curvature is large compared to all microscopic lengths 
in the problem so that it can be replaced by quasi-one dimensional wire along the 
x-axis with the twisted boundary conditions *&(L) = W(0) exp[2m(j>/<fio]. 

To solve this problem we need to take into account the ac perturbation only in 
the denominator q of the geometric series Ea.( l2.15t that determines the Cooperon. 
This is because 1 — q is governed by small corrections with the energy scale much 
smaller than 1 /r. The corrections due to the ac flux perturbation are given by the 
diagrams in Fig. 11. Note that in linear in V corrections one should take into 
account the small momentum k whereas in quadratic in V corrections k can be 
set zero. All the corrections of Fig. 1 1 can be computed by doing the pole integrals 
over £ and angular integrals over n as in the previous section. Note also that 1 — q 
is essentially the inverse Cooperon operator. For an arbitrary time dependence 
of ip(t) one should replace — iui in Eq. j2.18t by the time derivative d tl — d t . 2 
which one can obtain from the second term in Eq.( I3.12> . Such a structure of the 
time derivatives implies that the sum of time arguments t\ + t% = t[ + t' 2 is 
conserved [see Fig. 12] which is the consequence of the constant density of states 
approximation. As a result, the equation for a time-dependent Cooperon takes 
the form [5, 10, 11]: 

( 7T + ^ ^ + ^ + ^ - K] 2 } Ct(r,, rf; k x ) = 6{ * ~ f , (4.2) 

where the momentum \c x (cf)) — k m (<f>) — (27r/X)(m — 20/0o) is quantized 
according to the twisted boundary conditions. 

For completeness we give also an equation for the time-dependent diffuson 
[Fig.13]: 

jJ^+A) [<p(t + T]/2)-<p(t-v/2)-K} 2 \ D n (t,t';K) = 5 2^2 ><- 43) 
where k x = (2n/L) m is independent of the DC flux <p. 
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Fig. 1 1 . The ac flux and time-derivative corrections to the inverse Cooperon operator. 



In this case the structure of the time-derivative dt 1 + dt 2 suggests that the 
difference of the time arguments is conserved t\ — t% = t\ — t' 2 . 

In general the current response I(t) to the electric field E(t — r) is given by: 

I(t) = / K(t, t) E(t - t) dr, (4.4) 
Jo 

where in our particular case E(t) = —dt<p(t). 

At small 1/pf& the main contribution to the nonlinear response function is 
that of the quantum coherence correction shown in Figs. 5,7. All what we have to 
do in order to compute this nonlinear response function is to substitute the time- 
dependent Cooperon into Eq. j2.19l > and take care of all the time arguments [see 
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Fig. 12. The time arguments in the Cooperon. 
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Fig. 13. The time-dependent Diffuson. 



Fig. 14]. Namely, (i) according to Ea. d3.12t the time arguments corresponding 
to the beginning and the end of any solid line representing the disorder averaged 
functions {G^ A ) should be the same and (ii) the sum (difference)of "incoming 
times" in the Cooperon (Diffuson) are equal to the sum (difference)of "outgoing 
times" as in Figs. 12, 13. Note that application of rules (i),(ii) to the diagram of 
Fig. 14 leads to the coinciding time arguments t\ — > t 2 = t — t in the correspond- 
ing retarded-advanced junction (see Fig. 10). This ensures fulfillment of Ea. J4.4> . 
as ho(h - h) (V(h) - V(ti)) — v x dt<p\t - t) « E(t - t). 

The disorder average current response function K(t, r) is equal to: 

K(t, r) = a D r~ l e-^ - C t _ r/2 (r, -r), (4.5) 

where C t {r},rf) = {2t:vt 2 /Vol) X) kx CtOWlks) and Ct(i],f]';k x ) is the so- 
lution to Ea.d4~2l: 



2nvT 2 C t (ji,Ti';k x ) = 6(r, - rf) x 



(4.6) 
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V x d/dt (p (t -T ) 



Fig. 14. Nonlinear current response function and the Cooperon. 



exp 




/'[^ + C/2) + <^-C/2)-k x ] 2 d( 



} 



The first term in Eq.( l4.5> is just the linear response given by the diagrams Fig.lc,d. 
The second term corresponds to the quantum-coherent contribution of Fig. 14. 
The fact that C f _ T /2(r, — r) depends not only on r but also on t makes it possi- 
ble to have dc response caused by ac electric field E(t) = —d t <p{t). 

For a particular case of a ring with a constant and an ac magnetic flux we 
note that the dc response arises only from odd terms of expansion of Ea. (l4.6> in 
powers of the ac perturbation ip(t). As odd terms in ip(t) enter in a combination 
(p(t)\t x the dc response involves only odd terms in \s. x = k m = (2ir /L)(m — 
2(/)/<f)o), If the constant magnetic flux (/> = or we neglect quantization of mo- 
mentum and do an integral over k x instead of a sum, the result for the odd in 
ip(t) part of Ct{rj 1 77') is zero. Otherwise it is periodic in the flux <fi with the pe- 
riod 4>o/2, since fc m (</> + </>o/2) = fc. m _i(</>) and the summation is over all integer 
to. We see that the dc response to the ac perturbation, or the rectification of the 
ac flux by an ensemble of mesoscopic rings, is an essentially quantum, Bohm- 
Aharonov-like effect. Furthermore, it is an odd in the flux Bohm-Aharonov effect 
that can be represented by the Fourier series: 



Expression Eq.( I4.7> has the same symmetry and periodicity in the magnetic flux 
(j> as the disorder-averaged equilibrium persistent current [12, 13]. However, its 
magnitude may be much larger (in the grand-canonical ensemble considered here 
the disorder-averaged persistent current is strictly zero). 




(4.7) 
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Applying the Poisson summation trick 



(4.8) 



J dxf(x -(t>)^2< 



2-ninx 



dxe 2 * mx f(x) 



to Ea.( l4.6> one obtains 

AieD 



In = S ^T J drCl n) {T)d tV {t-T/2), 
where the overline means averaging over time t and 



IT 2 - p mSiM -tS 2 [v] 



(4.9) 



(4.10) 



Here td = L 2 /Do is time it takes for a diffusing particle to go around the ring, 
and Si, 2 are defined as follows: 



S 2 



2L 



t+r/2 



t-r/2 



2L(ip tl )t;r, 



(4.11) 



(4.12) 



Eqs.( l4.10l -( I4.12> are valid for an arbitrary time-dependence of <p(t). However, 
they take the simplest form for a noise-like ac flux with the small correlation 
time To <C td (but we assume To ^> r). In this case the second and the third 
terms in Eq.( l4.12> can be neglected and the first term reduces to a constant that 
determines the dephasing time caused by ac noise: 



— = 2D <pi(t). 



Then the time averaging in Eq. (l4.9> reduces to 



(4.13) 



id t <p(t- r/2)expi —2L 



t+r/2 



f(ti)dti 



2L(p 2 {t). (4.14) 



t-r/2 



Since the time-average of the a total time-derivative is zero, we can transfer the 
differentiation to the exponent. For the case of a white-noise ac flux only the 
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lower limit of integration in the exponent should be differentiated, the quantity 
(Vii )t;r being set zero afterwards. The remaining integral over t is done exactly: 



dr 

7m 



exp 



n 2 T D 
4r 



= 2 ^ e -«\/W^ 
riy/fB 



(4.15) 



I 




Fig. 15. The dependence of I\ on the ac power P 
perturbation. 



ip 2 (t) for a white-noise (a) and harmonic (b) 



Finally we arrive at a remarkably simple result for the disorder-average dc current 
generated in mesoscopic rings by a white-noise ac perturbation [8]: 



(4.16) 



4 






L ' 






exp 


—n— — 


7T 






Lip _ 



In 

where the dephasing length is equal to: 
1 1 



(4.17) 



It follows from Eq. j4.16> that the weak ac white noise produces a net dc current 
in an ensemble of mesoscopic rings which is of the order of e /t v where t v is 
the dephasing time produced by the same ac noise. At an ac power large enough 
to produce a dephasing length smaller than the circumference of a ring times 
the winding number n, the destructive effect of ac perturbation prevails and the 
current I n decreases exponentially [see Fig. 15]. 
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Note that the tail of the dependence /„ on the ac power is very sensitive to 
the correlations in the ac perturbation at different times. For instance, in the case 
of harmonic perturbation where ip(t) has an infinite range time-correlations, the 
current decreases very slowly, only as the inverse square-root of the ac power [6]. 
We will see later on that this is related with the phenomenon of no-dephasing 
points. 



5. Diffusion in the energy space. 

In this section we apply the Keldysh formalism outlined above to the problem of 
heating by external time-dependent perturbation. The main object to study will 
be the non-equilibrium electron energy distribution function: 

f(e;t) = l-~Jdr ] h t (r l )e- is o (5.1) 

which is related to the Keldysh function G K {t, t') = Vol^ 1 J d d r G K (i, t'; r, r) 
at coincident space variables (averaged over volume): 

G K {t + r l /2,t~n/2) = -2mvht{rj). (5.2) 

Eqs. (I5. 1I . J5.2I generalize Ea. J3.3> to the case of non-equilibrium, time-dependent 
energy distribution function. The total energy £ (t) of electron system can be ex- 
pressed in terms of f(e, t): 

S(t)=v Vol J dee [f{e, t) - 9(-e)} + const. (5.3) 

Then the time-dependent absorption rate W(t) = d t £(t) is given by: 

W(t) = lim dt d v G K = ? Km d t d v h t (n), (5.4) 

where 6 = (vVol)' 1 is the mean separation between electron energy levels 
(mean level spacing). 

For tutorial reasons in this section we consider a specific model system de- 
scribed by the Hamiltonian: 

H = e(p) + U(r) + V(r)<p(t), (5.5) 

where not only U(r) given by Eq. j2.1 H but also the perturbation potential V(r) 
is a Gaussian random field which is statistically independent of U(r) and is de- 
scribed by the correlation function: 

(V{r)V(r')) = — S(r - r'). (5.6) 
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This model is the simplest example of a potential ac source in contrast to Iheflux 
ac source considered in the previous section. In the low-frequency lot^ <C 1 
and the modestly low ac intensity Ft -C 1 limits this model is equivalent to the 
random matrix theory with the time-dependent Hamiltonian [5, 1 1]: 

H RM T=H + Vip{t), (5.7) 

where both Hq and V are random real-symmetric N x N matrices from the inde- 
pendent Gaussian Orthogonal Ensembles described by the correlation functions 

(CHo)„m) = (V nm ) = 0: 

({Ho)nm(Ho)n'm') = N(5/tt) 2 [5 nn >8 mm ' + 5 nm >5 mn >], (5.8) 

with S = \ j[v Vol) being the mean level spacing. 

In this RMT limit all the disordered and chaotic systems described by a real, 
spin-rorational invariant Hamiltonian are believed to have a universal behavior in 
the limit L, N — > oo which is characterized by only two parameters 6 and T and 
one dimensionless function ip(t). 

Let us first consider the disorder-averaged G K in the non-crossing approxima- 
tion. Since there are no vector vertices in the present model, the ladder diagrams 
(Lose Diffuson) similar to Fig. 2 make the main contribution. The corresponding 
diagrams are shown in Fig. 16 where the wavy line describes the (VV) correlator, 
Ea. J5.6l >. The condition 

Tt < 1 (5.9) 

allows to keep only the linear in V term in the expansion Ea. P.lOt for the 
(G R / A ) function attached to the RA-junction. Note that the diagrams Fig.l6a,b 
have opposite signs because J d£ (£ ± i/2r)~ 1 = ^firi. The result of calculation 
of both diagrams is: 

SG K = (27WT) 2 (i7w)(r/7rz/)Z^(i,i') (5.10) 
x {cp(t> + r 1 /2)-^(t' - V /2)} 2 h Q (r ] ). 

so that 

h t (r)) = (l-D r ,(t,t')T[cp(t' + r l /2)-Lp(t' -ri/2)} 2 ^) h (r,), (5.11) 

where the integration over t' is assumed; ho(r]) is determined by Eq.( l3.9l l. and 
D v (t, t') = 2-kvt 2 D v (t, t'; k = 0) is given by Eq.lgj) with D Q V. One can 
check using this equation that h t (rj) obeys the equation: 



{d t + T [<p{t + V2) - <p(t - r,/2)] 2 } h t (rj) = 0. 



(5.12) 
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t + T/2 



t -T/2 



R 



(p(f + T/2 ) 




t ' - T/2 



t + T/2 



R 



t -T/2 



A 



t' + T/2 




(p(t'-T/2 ) 



Fig. 16. The Loose Diffuson diagrams for the absorption rate: the dotted lines represent the (UU) 
correlation function while the wavy line corresponds to the { VV) correlator. 



This allows to give an explicit solution for h t (r]): 



h t (rj) = h (r]) exp 



b(C + W2)-^(C-W2)] 2 rfC 



(5.13) 



Note that a nontrivial dynamics of h t (rj) and thus that of f(e;t) is hidden in 
the exponential factor which is a purely classical object. The Fermi statistics of 
electrons considered here is taken into account by an initial distribution ho(r)). At 
t]uj <C 1 where w is the typical frequency of oscillations in (p(t), or equivalently 
at the energy resolution 5e ^S> u>, one can approximate (p(^+r]/2) — (p((~T]/2) « 

(dfip) rj. Then the exponent in Ea. (l5.13> reduces to exp —t T (dtip) 2 r] 2 
corresponds to Ea. d5.12> of the form: 



which 



{d t + D E ri 2 } h t (v) = 0, {d t -D E d 2 } f(e,t) = 0. 



(5.14) 
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In this approximation we come to the diffusion equation for the energy distribu- 
tion function with the energy diffusion coefficient: 



D E =T(d t <p) 2 ~Tuf. (5.15) 

For a harmonic perturbation ip(t) = cos(uit) the form of energy diffusion coef- 
ficient allows a simple interpretation. This is a random walk in the energy space 
due to a sequential absorption or emission of photons with the energy hw. Then 
the size of an elementary step is ±%uj and the rate of making steps is F/%. Now it 
is clear that the condition Ea. J5.9> that allows to use the approximation Ea. d3.10> 
has a physical meaning of a condition to absorb/emit at most one photon dur- 
ing the time of elastic mean free path. The opposite condition would mean that 
many photons can be absorbed/emitted during the time r which implies inelastic 
processes being stronger than the elastic ones. 

Using the diffusion equation we express dtht(rj) = — D e rj 2 h t (ri) . Then 
Ea. (l5.4> gives for the net absorption rate: 



W(t) = -in {D E /S) lim d^rfhtiv))- (5-16) 

■q— >0 

It is important that 

, r?->0. (5.17) 

TTT] 

for all ht(rj) corresponding to the Fermi-like energy distribution function f(e;t) < 
1 with f(e — > +oo; t) = 0, f(e — > — oo; t) = 1. Thus the existence of the Fermi 
sea makes the net absorption rate non-zero despite the diffusion character of the 
energy distribution dynamics: 

W(t) = W = -f-. (5.18) 



6. Quantum correction to absorption rate. 

We see that the non-crossing, or zero-loop, approximation leads to the classical 
picture of a time-independent absorption rate, the so called Ohmic absorption. 
This is in line with the fact established below that the same approximation leads 
to the classical Drude conductivity. Let us go beyond and consider the diagrams 
for G K that contain one crossing or one Cooperon loop. Different ways of pre- 
senting such a diagram are shown in Fig. 17. This is essentially the same fan 
diagram with the Loose Diffuson attached to it. However, because the perturba- 
tion does not contain a vector vertex, one should take a special care about the 
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Hikami box. It is shown in Fig. 18. Here a new element - besides the correlation 
between two V vertices that first appeared in Fig. 16- is the dotted line installed 
between two retarded or two advanced Green's functions. The rules of such an 
installation are (i) non-crossing of dotted-dotted and dotted-wavy lines and (ii) 
the presence of both retarded and advanced functions in each "cell" separated by 
dotted or wavy lines (otherwise the ^-integral zero). One can see that the sum of 
three diagrams in Fig.l8a-c and in Fig.l8d-f are zero. The remaining diagrams 
Fig.l8g and Fig.l8h have opposite signs so that the R-A junction and the V(t) 
vertices (denoted by an open circle in Fig. 18) that appear as the first term of 
expansion of (G R ^ A ) compose a combination 



T {<p(t' + ?7/2) - <p(t' - 77/2)) (<p(t' - h) - <p(t' -h- r;)) h Q ( V )(6.l) 

« rdMt')d t ,<p(t' -t^ifh^ri). 

The contribution to the absorption rate that comes from the diagram of Fig. 17 
is computed with the help of Ea. (l5.4> . The limit lim j; _>o d n r) 2 ho(ri) results in a 
universal constant i/w, and the time derivative dt of the Loose Diffuson gives a 
S(t — t') function in the leading order in V. So we obtain a quantum correction 
to the absorption rate [14]: 

' )U 1/1 - ' " f at^(t)^v(t-*i)C , t-t 1 /2(ti > -ti)«i- (6-2) 



W ttD e J 
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Fig. 18. Hikami box for a scalar vertex. 



where the limits of integration are fixed by an assumption that a time-dependent 
perturbation has been switched on at t = 0; Ct(rj, if] k) is the Cooperon, and 

Ctfarf) = (2^T 2 )^- l J2Ct(v,v'^)- (6.3) 

k 

Here we have to make some remarks. Ea. J6.2> is general and gives the first quan- 
tum correction to the absorption rate in a disordered or chaotic system of any 
geometry. However, the quantity C t (rj,ri') depends on the type of perturbation 
and on the system geometry. For a closed system of finite volume Vol there al- 
ways exists a regime where geometry of the system does not play any role. This 
is the limit where the typical frequency uj of perturbation is much smaller than the 
so called Thouless energy Exh- In disordered systems with the diffusion charac- 
ter of electron transport, the Thouless energy coincides with the gap between the 
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zero diffusion mode that corresponds to k = and the first mode of dimensional 
quantization that corresponds to k ~ 2ir/L , where L is the largest of the system 
sizes. In this case it is of the order of the inverse diffusion time Eth ~ 1/td. For 
such small frequencies one can neglect the higher modes with nonzero values of 
k and consider the zero diffusion mode with k = that always exists in closed 
systems. In this zero-mode, or ergodic limit the actual shape of the system does 
not matter at all. One can show that this is exactly the limit where the results 
obtained using the Hamiltonian Eq.( I5.5> are equivalent to the results obtained 
starting from the random matrix theory Ea. J5.7> . 

However, even in the ergodic limit the results for the Cooperon depend on the 
topology of the system. The equation for the Cooperon Ea. (l4.2> that corresponds 
to the non-single connected toplogy of a ring and a global vector-potential per- 
turbation differs from that for the scalar potential perturbations with local corre- 
lations in space described by Eas. d5.5l .( l5.6l . 

To see the difference we re-derive the Cooperon for the case of a scalar po- 
tential perturbation. The corresponding diagrams analogous to those shown in 
Fog. 1 1 for the case of the global vector-potential perturbation are given in Fig. 19. 
As a result of calculation of the corresponding ^-integrals we obtain an equation 
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Fig. 19. The ac scalar potential corrections to the inverse Cooperon. 



(in the ergodic limit k=0): 

1.2-^- + T [ip{t + V2) - <p(t - ?7/2)] 2 j C t ( V , rf) = 2% - rf). (6.4) 

This equation contains the difference tp(t + 77/2) — ip(t — T]/2) instead of the sum 
in Ea. (l4.2> . One can show [5] that the time-dependent random matrix theory that 
reproduces the Cooperon with the sum cp(t + rj/2) + ip(t — ry/2) corresponds 
to the random perturbation matrix V which is pure imaginary anf/-symmetric, 
rather than the real symmetric as in Eq. (I5.7> . 
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The corresponding equation for a time-dependent Diffuson appears to be the 
same as Ea. J4.3> : 

|^ + T [<p(t + 7?/2) - p(t - 77/2)] 2 | D v (t, t') = 5{t - t'). (6.5) 



To complete the solution Ea. J6.2> we give a solution to Ea.( l6.4> which is valid 
for an arbitrary time dependence of ip(t): 



r 



C t ( V , rf) = 0(V - V') exp { -- I [p{t + C/2) - <p(t - C/2)]' |> ■ (6.6) 



7. Weak dynamic localization and no-dephasing points. 

In this section we concentrate on the oscillating time dependence of perturbation 
with ip(t) = and (p 2 (t) = 1. The simplest example is that of a white-noise with 
ip(t)ip(t') = for f ^ t'. Then Eq.i6"6l takes the form: 

C t ( V , 77') = 0(77 - 77') exp [-T (77 - 77')] . (7.1) 

The negative exponential factor describes dephasing by a white noise perturba- 
tion. Because of the factor dt(fi(t)dt(fi(t — t\) the effective range of integration 
in Ea. J6.2t is of the order of the correlation time tq -c 1/T. This makes the 
correction to the absorption rate vanishing 6W(t)/Wo — * in the white-noise 
limit To — > 0. 

Now let us consider the simplest example of an infinite — range time corre- 
lations. This is the case of a harmonic perturbation ip(t) — cos(ojt). Then for 
U)\rj — rj'\ ^> 1 the Cooperon takes again the form Eq. (I7.1> but with the time- 
dependent dephasing rate: 

rf^rsin 2 ^ (7.2) 

One can see that there are certain moments of time uot n = vrn (n = integer) 
where the dephasing rate is zero. We will refer to these points t„ as no-dephasing 
points [15]. 

This remarkable phenomenon is generic to harmonic perturbation and is also 
present for the case of an ac vector-potential considered in the previous section. 
Its physical meaning is quite simple. Consider again a pair of electron trajectories 
with a closed loop shown in Fig. 8. In the presence of a time-dependent vector- 
potential A(t) the phase difference between these two trajectories is no longer 
zero but a random quantity: 

S$= [ A(t')[vi(0-v 2 (*')]^', (73) 
Jo 
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where v\^(t) are electron velocities at a time t and T is the time it takes for 
an electron to make a loop (the traversing time). For loops with the opposite 
directions of traversing vi(t') = — v 2 (T — t') [see Fig. 20]. Then we obtain 









' \(T-t') 



Fig. 20. Trajectories with loops traversed in opposite directions and synchronization. 



<S$ = / ' [A(T/2 + t')+A(T/2-t')]v 1 (t' + T/2)dt' (7.4) 

J-T/2 

Now assuming the period of oscillations of A(t) to be large compared to the 
velocity correlation time r we obtain after averaging over vi: 

((«5$) 2 } = D [ [A(T/2 + t')+A(T/2-t')} 2 dt'. (7.5) 

J-T/2 

We see that in the presence of a harmonic vector-potential A(t) oc sin cot the 
phase difference 6^ — vanishes for all the loop trajectories with the traversing 
time equal to the integer period of the time-dependent perturbation: 

27T71 

T= , n= 1,2,3... (7.6) 

u> 

Such loop trajectories with the traversing time synchronized with the period of 
perturbation play a crucial role in all the quantum coherence (interference) phe- 
nomena in disordered and chaotic systems (their effect on the universal conduc- 
tance fluctuations has been studied in Ref. [15]). Any periodic in time pertur- 
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bation leads to a selection of trajectories: for strong enough perturbation all tra- 
jectories that do not obey Ea. il. 6> do not interfere because of the random phase 
difference <5$ ^> 1. 

In the particular case of the quantum corrections to the energy absorption rate 
Ea. J6.2t the time-dependent Cooperon 

Ct- tl /2(h,-h) = exp [-2tir t _ tl/2 ] (7.7) 

is small everywhere except for the vicinity of no-dephasing points 

ujt[ n} = 2ut - 2m, n = 0,±l,±2... (7.8) 

For Ri ~ Tt > 1 one can expand r t _ tl/2 ~ {u/2) 2 T(t 1 - t^ f in the 
vicinity of no-dephasing points, do the Gaussian integration over C, = (ti — 4™') 
and put t\ = ti elsewhere. In particular d t if{t — t±) = —d t ip{t). Then we 
arrive at: 



SW(t) 

~wT 



-(2sin 2 ^)^/ d(e- D ^ 2ti "\ (7.9) 

^ _ J — oo 



Replacing ~^2 n by the integral (w/2n) J Q dt\ and averaging over time intervals 
much larger than 2tt/uj we finally get: 

Wit) Suj f f , f + °° dC n r 2 t IT 

= 1 / dh — — - e _ ^ = 1 — \ — , (7.10) 



(7.11) 



W tt J J_„ 2tt V t 

where 

** ~ 26* ■ 



We have obtained a remarkable result that the absolute value of the quantum 
correction to the absorption rate is growing with time. This is the consequence 
of the existence of no-dephasing points, as otherwise the exponentially decaying 
Cooperon would lead to a saturation of the integral over t\ in Eq. j6.2t at large 
times t. The negative sign of the correction implies that the absorption rate, or 
energy diffusion coefficient decreases with time. This phenomenon can be called 
"weak dynamic localization" in full analogy with the weak Anderson localization 
when the diffusion coefficient in space is decreasing with the system size. It 
has been first discovered [16-18] for a simple quantum system - quantum rotor 
subject to the periodic (5-function perturbation. At t ~ t* the quantum correction 
is of the order of the classical Ohmic absorption, and we can expect the strong 
dynamic localization to occur. 
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Eq.( I7.10> suggests a more precise relationship between the dynamic localiza- 
tion for a quantum system in the ergodic (zero-dimensional) limit subject to a 
harmonic perturbation and the Anderson localization in a quasi-one dimensional 
disordered wire. We observe that 

J-oo 2tt J_ qo 2vr J_ QO 2vr D E ( 2 ~ iuj' 

At the same time, according to Ea. J2.22U he correction the (complex) frequency 
dependent conductivity of an infinite quasi-one dimensional wire is proportional 
to: 

, , . f + °° dk 1 

<H*0 = / ^ n , 2 ■ ■ (7-13) 

We see that the deviation from the no-dephasing points £ = t\ — is an ana- 
logue of the momentum k. Then with a suitable choice of a time and frequency 
scale one obtains: 

W 2vr (-iw' + 0) \ a )' ' 

Note that Eq.( I7.14> is non-trivial as it establishes a relationship between an essen- 
tially non-equilibrium property of a zero-dimensional system and an equilibrium 
property of a quasi-one dimensional system. 

This relationship can be proven for any diagram with an arbitrary number of 
the Diffuson-Cooperon loops (for the case of two loops see Ref. [19]) and we 
conjecture that it is exact in the ergodic (or random matrix) regime for tT, tui 3> 
1. 

Eq.( I7.14> helps to establish a character of decay of the absorption rate W(t) 
for t ^> i*. To this end we recall the Mott-Bereznskii formula for the frequency- 
dependent conductivity in the localized regime [20]: 

cr(w) cx uj 2 In 2 uj. (7.15) 

Then using Ea. d7.14> one obtains the corresponding absorption rate W(t) &tt^> 
i ~ * • 

W(t) cx ^. (7.16) 
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W(t)/W f 




Fig. 21. The time dependence of the absorption rate. 

8. Conclusion and open questions 

The goal of this course was to give a unified picture and a unified theoretical tool 
to consider different quantum coherence effects in disordered metals. The unified 
picture is that of interference of electron trajectories with loops traversed in the 
opposite directions. The corresponding theoretical machinery is the Diffuson- 
Cooperon diagrammatic technique. We have demonstrated that even the first 
quantum correction diagram with the Cooperon loop can describe such different 
and nontrivial phenomena as weak Anderson localization, quantum rectification 
and dynamic localization in the energy space. We did not try to give a review of 
the development in the field but rather to concentrate on few important examples 
and to demonstrate how does the machinery work in different cases. That is why 
many related issues have not been discussed and the corresponding works have 
not been cited properly. We apologize for that. 

There are few open problems that are related with the main subjects of this 
school and in our opinion are warrant a study. This is first of all a unified the- 
ory of energy absorption where both the Zener transitions picture [21] and the 
sequential photon absorption picture are incorporated. The suitable theoretical 
tool for that is believed to be a nonlinear sigma-model in the Keldysh repre- 
sentation derived in Refs. [14]. The perturbative treatment of this field theory 
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reproduces the diagrammatic technique discussed in this course, and the non- 
perturbative consideration in the region r, u> < 5 should give the results obtained 
in the framework of the Zener transitions picture. 

Another possible direction is the role of interaction for the quantum rectifica- 
tion and dynamic localization. Some of the interaction effects in dynamic local- 
ization have been recently considered in Refs. [22-24] using the Fermi Golden 
Rule approximation. However, an interesting regime of localization in the Fock 
space [25] where the Fermi Golden Rule does not apply is awaiting an investiga- 
tion. 
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